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The c(2x2) reconstruction of (001) PbTiC>3 surfaces is studied by means of first principles calcu- 
lations for paraelectric (non-polar) and ferroelectric ([001] polarized) films. Analysis of the atomic 
displacements in the near-surface region shows how the surface modifies the antiferrodistortive 
(AFD) instability and its interaction with ferroelectric (FE) distortions. The effect of the surface 
is found to be termination dependent. The AFD instability is suppressed at the Ti02 termination 
while it is strongly enhanced, relative to the bulk, at the PbO termination resulting in a c(2x2) 
surface reconstruction which is in excellent agreement with experiments. We find that, in contrast 
to bulk PbTiOs, in- plane ferroelectricity at the PbO termination does not suppress the AFD in- 
stability. The AFD and the in-plane FE distortions are instead concurrently enhanced at the PbO 
termination. This leads to a novel surface phase with coexisting FE and AFD distortions which is 
not found in PbTiQ 3 bulk. 



PACS numbers: 

I. INTRODUCTION 

The ABO3 perovskite structure exhibits several lat- 
tice instabilities which lead to a rich variety of phase 
diagrams with structures ranging from non-polar anti- 
ferrodistortive (AFD) to ferroelectric (FE) and antifer- 
roelectric (AFE). In these systems, polar and non-polar 
instabilities often compete and tend to suppress one an- 
other. For instance, in SrTiOa it has been found that 
ferroelectricity is enhanced if the AFD degrees of free- 
dom are artificially frozen out by setting the correspond- 
ing distortion to zeroi^ Of these various instabilities, 
the polar instabilities are the ones responsible for the 
technologically important and fundamentally interesting 
properties, such as ferroelectricity and large piezoelectric 
and dielectric responses. 

Currently there is a great deal of effort being expended 
to develop nanoscale devices. If nanoscale systems based 
on ferroelectricity are to be feasible, then the polar in- 
stabilities must remain dominant down to the nanoscale. 
In systems of reduced dimensionality, such as thin films 
or nanoparticles, a substantial fraction of the material is 
found in proximity to a surface, and surface effects may 
become important. The surface can affect the material 
behavior by modifying the strength of the various insta- 
bilities and their interactions leading to novel phases not 
present in bulk systems. 

Lead titanate (PbTiOs) is the prototype of a large class 
of Pb-based perovskites. The cubic perovskite struc- 
ture of PbTiOa possesses several branches of unstable 
phonons, including FE instabilities which consist of zone- 
center polar TO modes and AFD instabilities which con- 
sist of nonpolar zone-boundary modes involving rotations 
of oxygen octahedral cages surrounding the Ti atoms^ In 
bulk PbTi0 3 , the FE and AFD instabilities of the cubic 
structure compete with each other and the FE lattice 
distortion suppresses the AFD distortion so that AFD 
distortions do not participate in the ground-state tetrag- 



onal FE structure. The possible impact of the proximity 
of the surface on lattice instabilities and their mutual 
interactions is thus of great importance in understand- 
ing the structure and properties of PbTiOa surfaces and 
ultrathin films. 

Previous theoretical studies on lead titanate have been 
confined to the behavior of the polar FE instability at the 
(001) surface. 4 ' 5 ' 6 These studies employed (lxl) surface 
periodicity, which by construction does not permit AFD 
deformations. Proximity to the surface, however, could 
enhance the strength of non-polar instabilities and mod- 
ify their interaction with the FE instability. Indeed, an 
AFD c(2x2) reconstruction has recently been found at 
the (001) surface of PbTiOa using grazing incidence x- 
ray scattering^ suggesting that the AFD instability is 
important at the surface. 

In the following, we study the behavior of the AFD 
instability and its possible competition with ferroelec- 
tricity at the PbTiOa (001) surface using first-principles 
calculations. Initially we study the AFD distortion at the 
surface, using symmetry constraints to freeze out polar 
distortions; next we relax this constraint and study the 
competition between the AFD and in-plane FE distor- 
tions at the PbO-terminated surface, which is the equi- 
librium termination for PbTiOs (001) 5 ". 

In the absence of ferroelectricity, we find that the im- 
pact of the surface is termination-dependent: at the PbO- 
terminated surface, the strength of the AFD instability 
is enhanced relative to the bulk; while at the Ti02 ter- 
minated surface, the strength of the instability is dimin- 
ished. At the PbO-surface, in contrast to the bulk, where 
the presence of ferroelectricity suppresses AFD, we find 
that the two coexist and we observe a c(2x2) reconstruc- 
tion that agrees very well with experiment^ 

The paper is organized as follows: section[n]is devoted 
to the theoretical method and the technical details of our 
calculations. In section ILTTI we study the AFD instability 
in bulk cubic lead titanate. In section ITVl we study the 
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surface effect upon the AFD distortion in absence of FE 
distortion. The competition of AFD and in-plane FE dis- 
tortions in the near-surface region is studied in sectionlvl 
The conclusions are summarized in section Wl\ 

II. THEORETICAL METHOD 

Our theoretical study is based on ab initio calcula- 
tions that were performed within density-functional the- 
ory (DFT) using the plane-wave pseudopotential method 
as implemented in the PWSCF package^. In this section 
we describe the technical details of the calculations and 
the supercells used to study the surface of PbTiC>3 (001) 
films. 



A. Computational details 

The exchange and correlation energy was given within 
the local density approximation (LDA) using the pa- 
rameterization of Perdew and Zunger 9 , Ultrasoft 
pseudopotentialsiS were used to describe the electron-ion 
interaction, treating as valence states the 3s, 3p, 3d, and 
4s states of Ti, the 4s, 4p, Ad, and 5s states of Zr, the 5d, 
6s, and 6p states of Pb, and the 2s and 2p states of O. The 
wave- functions were expanded in plane waves with an en- 
ergy cutoff of 25 Ry. A larger cutoff of 250 Ry was used 
for the charge density, which includes the augmentation 
charges required by the use of ultrasoft pseudopotentials. 
The Brillouin zone (BZ) integration was performed using 
a (3,3,2) Monkhorst-Pack meshii in the Brillouin zone of 
the c(2x2) supercells. Relaxations of atomic coordinates 
have been iterated until the forces on the atoms were less 
than 0.001 Ry/a.u. (26 meVA" 1 ). 

B. Surface and slab geometries 

In PbTiC>3, the stacking along [001] consists of al- 
ternating TiC>2 and PbO atomic layers (see Fig. [IJ. 
There are two possible terminations for the non-polar 
(001) surface: the PbO-terminated surface and the TiC>2- 
terminated surface. We have studied both types of sur- 
face termination. 

The surfaces were modeled using a periodic slab geom- 
etry. Specifically, we used symmetric slabs, terminated 
by two equivalent surfaces, and separated by a vacuum 
region of approximately 14 A. For the PbO-termination 
we used slabs of 11 atomic layers (7-layer slabs were used 
for direct comparison with previous calculations), and for 
the TiC>2 termination we used slabs of 9 atomic layers. A 
c(2x2) in-plane supercell was used, containing 4 atoms in 
each PbO layer and 6 atoms in each Ti02 layer (Fig.^l. 
The in-plane lattice parameter is a>oV2, where ao=7.37 
a.u. is the theoretical lattice parameter of cubic PbTiOa 
bulk. This mimics the epitaxial constraint imposed by a 
lattice-matched cubic substrate. 



In order to investigate the behavior of the AFD distor- 
tion in the absence of FE distortions, we studied non- 
polar (paraelectric) c(2x2) films. Specifically, we re- 
laxed the atomic positions in the lowest symmetry space 
group containing inversion (P2/m)^£ as this allows for 
the largest number of degrees of freedom for the atomic 
displacements during the energy minimization while con- 
straining the spontaneous polarization to zero. With this 
symmetry, all the atoms are allowed to move in the z (sur- 
face normal) direction and the oxygen atoms in the Ti02 
layers are also allowed to move in the xy plane (parallel 
to the surface). All the other atoms are fixed by symme- 
try to the ideal perovskite coordinates in the xy plane. 
All the symmetry-allowed displacements were fully re- 
laxed to find the equilibrium structure which minimizes 
the energy. 

The competition between the AFD and FE distortions 
was studied using a ferroelectric film with in-plane po- 
larization. This is achieved by using a symmetric c(2x2) 
slab with space group Pmii 3 . The mirror symmetry, M z , 
prevents spontaneous polarization in the z direction and 
allows only an in-plane polarization to develop in the 
film. 

The case of perpendicular polarization is not consid- 
ered in this work because it requires special treatment 
of the electrical boundary conditions and a separate 
discussioniia 




FIG. 1: (a) Crystal structure of the ideal cubic perovskite 
unit cell of PbTiOa. In (b) and (c) are top views of the ideal 
PbO and Ti02 (001) atomic layers; the solid lines indicate 
the c(2x2) surface unit cells. 

The atomic structure was optimized by relaxing the 
symmetry-allowed atomic displacements. We defined 
atomic displacements, S a , relative to the ideal perovskite 
structure. Atomic relaxations normal to the surface, 5 Z , 
cause a rumpling of the atomic layers and a change in 
interlayer distances relative to the ideal cubic structure. 
The in-plane atomic relaxations, <5||, chiefly consist of ro- 
tations of the oxygen squares surrounding the Ti atoms 
in the Ti02 planes. 

To analyze the optimized structures and compare with 
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previous results we computed the following parameters: 
the rumpling amplitude, 77 j , in each layer, i, the change in 
the interlayer distances relative to the ideal cubic struc- 
ture, Adij, and the rotation angle of the oxygen squares 
in the TiC>2 planes, 6i. 

The structural parameters are defined as follows. For 
each atomic plane, i, we let S z (Oi) and 5 z (Mi) be the 
average z displacements of the oxygen and metal atoms, 
respectively. As in Ref. |5j, the interlayer distance dij is 
given by the difference of the averaged S z displacements, 
[§z(Oi) + S z (Mi)]/2, computed for layer i and j. The 
layer rumpling is given by % = \5 z (Mi) - S z (Oi)}. A 
negative rumpling, 77, < 0, means that on average the 
metal atoms are deeper inside the surface relative to the 
O atoms. The rotation angles 0, are zero for the (lxl) 
surface, a value different from zero indicating a c(2x2) 
reconstruction. 



III. ANTIFERRODISTORTIVE INSTABILITY 
IN BULK CUBIC PbTiO s 

To understand how the AFD instability is modified 
near the surface we first study the AFD distortion in 
bulk PbTiC>3. The rotation of oxygen octahedra is a 
common AFD instability in perovskite oxides. In bulk 
cubic PbTiC>3 it is associated with two unstable zone- 
boundary phonons: the R25 mode at the R point of the 
Brillouin zone, q=(lll)7r/a, and the M3 mode at the M 
point, q=(110)7r/a. The three-fold degenerate R25 mode 
can be chosen to correspond to the rotations, about the 
(001) Cartesian axes, of TiC>6 octahedra in opposite di- 
rections from one cubic unit cell to the next. The singly- 
degenerate M3 mode is similar to the R25 mode, but 
neighboring planes of octahedra along the rotation axis 
rotate in-phase instead of out-of-phase. 

We considered distorted structures obtained by sepa- 
rately "freezing in" R25 and M3 modes in cubic PbTiC>3 
at the theoretical cubic lattice constant (ao=7.37 a.u.). 
In Fig. [3 we show the computed values of the total en- 
ergy per 5-atom unit cell as a function of the rotation 
angle. The energy is lowered by the AFD distortion, 
as expected, and the R25 distortion with rotation along 
a Cartesian axis is more favorable than the related M3 
distortion. The computed equilibrium value of the octa- 
hedral rotation angle is #buik = 3.3°, and the energy gain 
associated with the AFD instability in the bulk is AE^pp 
= 1.2 meV per bulk unit cell (well depth in Fig. [5J. 

The small energy gain indicates that cubic PbTiOs, in 
the absence of any FE distortion, displays only a weak 
AFD instability. In particular, it is much weaker than 
the ferroelectric instability, which has a well depth of 
about 15 meV for the cubic undistorted unit cell (and 
about 50 meV per unit cell if coupling with the tetragonal 
distortion is included). 

The AFD instability is suppressed by the FE distor- 
tion. In fact the ground-state ferroelectric tetragonal 
structure with 5 atoms per unit cell has no unstable 



phonons at any wavevectoni 5 * The ferroelectric structure 
results from the coupling of a FE polar distortion (TO 
mode) and a tetragonal strain. We checked, therefore, 
if the suppression of the AFD instability can be caused 
by the FE distortion alone. To do this we first relaxed 
the ferroelectric distortion without coupling to tetrago- 
nal strain by keeping the unit cell fixed to the ideal cu- 
bic structure and allowing internal atomic displacements 
along the (001) direction. As a result a spontaneous po- 
larization develops along the (001) direction. We then 
computed the phonons for this "cubic polarized" struc- 
ture (C-FE) and found no unstable modes at the R point 
of the Brillouin zone. We conclude that the FE distor- 
tion suppresses the AFD instability even in absence of 
tetragonal strain. 




0° 1° 2° 3° 4° 5° 6° 7° 8° 
oxygen-octahedra rotation angle 

FIG. 2: Energetics of the AFD instability in bulk cubic 
PbTiC>3. Filled diamond indicates the minimum, 9 = 3.3° 
and AE = 1.2 meV per bulk unit cell, to be compared with 
the ferroelectric well depth, AE = 50 meV per bulk unit 
cell. Triangles and circles are calculations for the M3 and 
R25 distortions, respectively. The solid line is a polynomial 
interpolation. 



IV. PbTiO-3 (001) SURFACES IN THE ABSENCE 
OF FERROELECTRIC DISTORTIONS 

To understand the behavior of non-polar distortions 
near the surface, we study the atomic structure of the 
non-polar (paraelectric) films. We studied both the PbO 
and Ti02 terminations with c(2x2) in-plane transla- 
tional symmetry. Our primary focus in this section will 
be on characterizing the c(2x2) reconstruction and its 
dependence on termination. 

A. PbO termination 

Results for the atomic displacement and structural 
parameters of the optimized c(2x2) surface structure, 
shown in Fig. [21 are presented and compared with those 
for the (lxl) surface in Tables HI and ITT1 

The results for the (lxl) surface are in excellent agree- 
ment with a previous calculation 5 ", given in Table [H] 
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Small differences are chiefly due to the different slab 
thickness used in the two calculations (we obtained al- 
most identical results using a 7-layer slab as in Ref|3). 
We used a thicker 11-layer slab to study the behavior of 
lattice distortions in the near-surface region and their de- 
cay to bulk-like behavior in the central layers of the slab. 
The change in interlayer distance displays an oscillating 
behavior: the first interlayer distance, d\2, is substan- 
tially contracted, whereas the second, c?23, is expanded 
and the third, d 3 4, is again contracted. The most no- 
ticeable features is the large inward displacement of the 
surface Pb atom relative to the surface O atom. This 
accounts for the large first interlayer contraction, Adi2, 
and first layer rumpling, r/i. 

TABLE I: Atomic displacements, relative to the ideal cubic 
perovskite structure, for the PbO-terminated (001) surface of 
non-polar PbTiC>3 films. The relaxations perpendicular (8 Z ) 
and parallel (5m) to the surface are given as a percentage of 
the cubic lattice parameter ao. 



layer 


atom 


c(2x2) 

S z 




(lxl) 

8 Z 


1 


Pb 


-1.83 




-4.34 







-0.39 




-0.41 


2 


Ti 


4-2.79 




4-2.67 







+1.85 


10.03 


4-1.26 


3 


Pb 


-1.57 




-2.08 







4-0.41 




-0.35 


4 


Ti 


+0.51 




4-0.64 







4-0.09 


2.52 


4-0.15 


5 


Pb 


-0.25 




-0.56 







-0.02 




-0.11 


6 


Ti 


















3.43 






The optimized atomic structure for the PbO termina- 
tion of the non-polar c(2x2) PbTiC>3 (001) film is shown 
in Fig. |3 displaying the large c(2x2) reconstruction oc- 
curring in the Ti02 subsurface layer. The reconstruction 
consists of alternating clockwise and anti-clockwise rota- 
tions of the Ti-centered oxygen squares in the TiC>2 layer. 
The in-plane oxygen displacements are 10% of the bulk 
lattice parameter, corresponding to a rotation angle of 
#1=11°. This atomic distortion is essentially an AFD ro- 
tation of the surface oxygen-octahedra about the z axis; 
particularly noteworthy is the substantial enhancement 
of the AFD distortion at the PbO-terminated surface. 
Neighboring planes of octahedra along the rotation axis 
rotate out of phase, as in the dominant AFD bulk insta- 
bility (R.25 mode). Deep inside the film, the AFD dis- 
tortion recovers its bulk value in the absence of a ferro- 
electric distortion. The relaxation energy, relative to the 
ideal (lxl) surface, is —0.29 and —0.18 eV per (lxl) sur- 
face unit cell for the c(2x2) and (lxl) surfaces, respec- 
tively. The c(2x2) reconstructed surface is much more 



TABLE II: Structural parameters for the PbO-terminated 
(001) surface of non-polar PbTiOs films. Change in the in- 
terlayer distance, Ady, and layer rumpling, rji, in percent of 
the lattice constant atmik- Rotation angle, 6i, of the oxygen 
squares in the Ti0 2 layers. Relaxation energy relative to the 
ideal (lxl) surface, AE rc i ax , in eV per (lxl) surface unit cell. 





this work 


Meyer et al. a 




c(2x2) 


(lxl) 


(lxl) 


Ad 12 


-3.4 


-4.3 


-4.2 


Ad 2 3 


4-2.9 


4-3.2 


+2.6 


Ad 3 4 


-0.9 


-1.6 


-0.8 


Ad 45 


4-0.4 


4-0.7 




Ad 5 e 


-0.1 


-0.3 




rji 


-1.4 


-3.9 


-3.9 


rp 


4-0.9 


4-1.4 


+ 1.2 


T]3 


-2.0 


-1.7 


— 1.2 




4-0.4 


4-0.5 




r\h 


-0.2 


-0.4 




02 


11.4° 










-2.9° 








e 6 


3.9° 







(#bulk) 


(3.3°) 






AE rc lax 


-0.29 


-0.18 




a Rcf 5 



favorable than the relaxed (lxl) surface, the difference 
in energy being 0.11 eV per (lxl) surface unit cell. It 
is interesting to notice that the symmetry of the relaxed 
structure (space group P4/nmm) is higher than the min- 
imal symmetry constraint imposed during energy mini- 
mization (space group P2/m). The four oxygen atoms 
in each TiC>2 layer are equivalent by symmetry in the 
relaxed structure. 

The main differences between the c(2x2) and (lxl) 
structures are (i) the AFD rotation of the Ti04-squares, 
chiefly in the second layer, (ii) a much smaller inward 
displacement of the surface Pb atoms, 6 z (Pbi) = —1.83 
% of ao, (hi) and, consequently, a substantial reduction 
of the first layer rumpling (771) and the first interlayer 
contraction (Adi2). 

These results can be understood in terms of the par- 
ticular chemistry of lead. In the ideal cubic perovskite 
bulk, Pb is in the center of an O12 cage and the twelve 
equivalent Pb-0 bonds are 2.76 A long. The tendency of 
lead to move off-center and form stronger covalent PbO 
bonds, as short as 2.59A, has been demonstrated to be an 
important factor in ferroelectricity in PbTi03ii§*i£ For- 
mation of shorter bonds is favored in the broken symme- 
try environment at the surface and, as we discuss next, 
can be regarded as the driving force for the c(2x2) AFD 
reconstruction occurring at the PbO-terminated surface. 



5 



a) 



b) 




FIG. 3: Optimized atomic structure of the c(2x2) non-polar 
(001) PbTi0 3 film with PbO surface termination, (a) Top 
view; only the first top two layers are shown, (b) Side view; 
only the top half of the 11-layer slab is shown. The oxygen 
atoms in the second layer are shown in a darker shade of gray. 
The c(2x2) AFD reconstruction, consisting of rotation of the 
Ti-centered oxygen squares in the Ti02 subsurface layer, is 
clearly visible. Only Pb-O bonds shorter than 2.6 A are shown 
(in black); they are 2.39 A long. 



In Table II I II we show the values of the interlayer 
bondlengths computed for the c(2x2) and (lxl) sur- 
faces. Each Pb atom forms four bonds with the in-plane 
oxygen atoms, four interlayer bonds with O atoms in the 
plane above, and four interlayer bonds with O atoms in 
the plane below; the surface Pb atom is undercoordi- 
nated and has eight bonds (four in-plane bonds and four 
interlayer bonds with the plane below). In the non-polar 
slab, in-plane Pb off-centering is not allowed by sym- 
metry. Thus, the four equivalent in-plane Pb-0 bonds 
can only change by bending; their lengths are essentially 
fixed to 2.76 A and are not included in Table ITTTI In con- 
trast, the interlayer Pb-0 bonds can stretch or contract 
substantially. When a (lxl) periodicity is enforced, the 
only possibility to shorten the Pb-0 bond, as found in 
Ref.||J is by a large downward displacement of the surface 
Pb atom, resulting in bondlengths (2.61 A long) com- 
parable to those in tetragonal FE bulk. At the c(2x2) 
surface the AFD distortion offers an alternative way, ev- 
idently more energetically favorable, to form fewer but 
shorter Pb-0 bonds. Each surface Pb atom has four 
oxygen neighbors in the subsurface Ti02 layer. Rota- 
tion of the Ti04-squares brings two of these much closer, 
resulting in even stronger bonds (2.39 A long) only a 
bit longer than those in bulk PbO (2.32 A long, experi- 
mental valued) . The rotation also shortens significantly 
two of the bondlengths for the Pb3 atom to 2.60 A , 
comparable to those in tetragonal FE bulk. Since the 
rotation provides a more effective alternative mechanism 
for shortening Pb-0 bonds, the inward displacement of 



surface Pb atoms (and the related rumpling rjx and in- 
terlayer contraction Adi?) is substantially reduced with 
respect to the (lxl) surface. The average bond length for 
each Pb atom does not change much between the c(2x2) 
and (lxl) surfaces and it is not much different from the 
undistorted bond length. 

As a further verification of the particular role of lead in 
the enhancement of the AFD distortion at the surface, we 
have relaxed the non-polar (001) surface of BaTiOa. The 
Ba-0 bond has a more ionic character and is therefore 
a good contrasting example. We find no AFD distor- 
tion at the BaO surface. This confirms that the forma- 
tion of stronger covalent bonds is the origin of the AFD 
reconstruction at the PbO-tcrminated (001) surface of 
PbTi0 3 . 

Finally we notice that, contrary to the Pb-0 bonds, 
the interlayer Ti-0 bonds are almost the same at the 
c(2x2) and (lxl) surfaces, and their values reflect the 
alternating contraction and expansion of the interlayer 
distances. 



TABLE III: Computed values for the Pb-0 and Ti-0 inter- 
layer bond lengths (in A ) at the PbO terminated c(2x2) and 
(lxl) surfaces of the non-polar (001) PbTi03 film. For com- 
parison, the in-plane Pb-O bond length is 2.76 A. Average 
bond lengths Pbi-O for each Pbi atom are given in the "aver- 
age" columns. The labels i and j indicate the atomic planes, 
with i,j=l being the first surface layer. Shorter Pb-O bonds 
are highlighted in bold. 







c(2 


x2) surface 


(lx 


1) surface 


Pbi 




short 


long average 




average 


Pbi 


-0 2 


2.39 


2.96 2.72 


2.61 


2.69 


Pb 3 
Pb 3 


-o 2 
-o 4 


2.60 

2.64 


3.13 2 7g 
2.78 


2.85 
2.70 


2.77 


Pbs 
Pb 5 


-o 4 
-o 6 


2.70 
2.66 


2.84 2 
2.85 


2.78 
2.74 


2.76 


Tii- 


-Oj 










Ti 2 


-Oi 


1.83 




1.83 




Ti 2 


-o 3 


2.04 




2.07 




Ti 4 


-o 3 


1.95 




1.91 




Ti 4 


-O b 


1.97 




1.98 




Ti 6 


-Ob 


1.95 




1.95 








B. 


Ti0 2 termination 







The computed structural parameters for the Ti02 ter- 
mination of the non-polar film are summarized in Ta- 
ble IIVI As we discuss in detail below, the large AFD 
surface reconstruction observed at the PbO termination 
does not occur at the Ti02 termination. 

Our results for the (lxl) Ti02 surface are in very good 
agreement with previous calculations. The structure of 
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the optimized c(2x2) supercell is little changed from that 
obtained for the (lxl) slab. The only difference is a small 
AFD rotation of the oxygen octahedra in the inner lay- 
ers, which is due to the AFD instability of bulk PbTiOs 
in absence of FE distortion. The AFD rotation in the 
first surface layer is essentially zero (6± = 0.1°) indicat- 
ing that, in contrast to the PbO termination, the AFD 
instability is weakened at the TiC>2 termination. The op- 
timized c(2x2) and (lxl) structures have the same en- 
ergy within numerical accuracy (they differ by less than 
10~ 6 Ry per atom). 

Mayer et al£ have shown that in thermodynamic equi- 
librium , the (lxl) PbO-termination is always more sta- 
ble than the (lxl) Ti02-tcrmination. Our results show 
that the Ti02-termination does not reconstruct while the 
PbO-termination is further stabilized by the AFD c(2x2) 
reconstruction. Therefore, the PbO termination remains 
the most favorable. 

The computed values for the Pb-0 and Ti-0 interlayer 
bond lengths are shown in Table [V] The short and long 
Pb-0 bonds are very similar due to the very small AFD 
rotation angles, and they are essentially equal at the sur- 
face where the rotation angle is zero within the accuracy 
of the calculation. The average bond length for each 
Pb atom is equal to the Pb-0 bond length in the ideal 
perovskite structure (2.76 A ). The interlayer Ti-0 bond 
lengths are very similar to those at the c(2x2) and (lxl) 
PbO terminations (see Table lrTT|l indicating that they are 
independent of both termination and reconstruction. 



TABLE IV: Structural parameters for the Ti02-terminatcd 
(001) surface of paraelectric PbTiOa films. Change in the 
interlayer distance, Adij, and layer rumpling, r]i, in percent 
of the lattice constant atmik- Rotation angle, 9i, of the oxygen 
squares in the Ti02 layers. Relaxation energy relative to the 
ideal (lxl) surface, AE rc i ax , in eV per (lxl) surface unit cell. 





this work 




Meyer et al° 




c(2x2) 


(lxl) 


(lxl) 


Adi2 


-4.3 


-4.3 


-4.4 


Ad 23 


+3.4 


+3.3 


+3.1 


Ad 34 


-1.2 


-1.3 


-0.6 


Ad 45 


+0.6 


+0.5 




Tjl 


-3.2 


-3.2 


-3.1 


m 


+4.3 


+4.3 


+4.1 




— l.U 


n q 
— u.y 


— u. * 




+0.9 


+0.9 




Is 


0.0 


0.0 




0i 


-0.1° 








03 


-2.1° 








05 


2.8° 







(0bulk) 


(3.3°) 
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C. Surface effect on the strength of the AFD 
instability 

The energy gain associated with the lattice deforma- 
tion is a measure of the strength of the instability, and 
can be used to compare the strength of the AFD insta- 
bility at the PbO-surface with the one in the bulk. 

The energy gain associated with the AFD deformation 
in a slab can be defined as the difference between the 
energy of the relaxed (lxl) non-polar slab, in which by 
symmetry there is no octahedron rotation, and the en- 
ergy of the relaxed c(2x2) non-polar slab. This energy 
difference can be divided into a bulk and a surface con- 
tribution as follows: 

2 £(lxl)slab _ £c(2x2)slab = jVA^ulk + 4A£^ (1) 

The bulk contribution, NAE%$^, is the energy gam as- 
sociated with the bulk-like rotation of the inner "non- 
surface" octahedra, with N being the number of such 
bulk- like octahedra in the c(2x2) unit cell (N=6 for the 
11 layer slab). The surface contribution, 4AE^p^, is the 
energy gain associated with the rotation of the surface oc- 
tahedra. The factor of four arises because in the unit cell 
of a c(2x2) slab there are four surface octahedra, two at 
each surface of the slab. Using Eq. 2] and the computed 
values for the three energies: £;( lxl ) slab ; £c(2x2)siab^ anc j 



TABLE V: Pb-O and Ti-O interlayer bond length in A com- 
puted for the Ti0 2 termination of the c(2x2) PbTi0 3 (001) 
surface. Average Pb-O bond lengths are also reported. We 
use the same notation as in Table Mil The in-plane Pb-O 
bond lengths are 2.76 A long. 



Pb, 


-Oi 


short 


long 


average 




Th- 




Pb 2 


-Oi 


2.62 


2.63 


2.76 


Tii 


-o 2 


1.80 


Pb 2 


-O a 


2.85 


2.95 




Ti 3 


-o 2 


2.02 


Pb 4 


-Oa 


2.67 


2.78 


2.76 


Ti 3 


-o 4 


1.90 


Pb 4 


-O s 


2.72 


2.86 




Ti 5 


-o 4 


1.95 



A£^pp , we find that the energy associated with the sur- 
face AFD distortion at the PbO-termination is AEf^= 
109 meV per (lxl) surface unit cell. 

By comparing this result with the AFD and FE well 
depths in the bulk (1.2 meV and 50 meV, respectively) we 
find that the energy gain associated with the surface AFD 
rotation is two orders of magnitude larger at the PbO- 
tcrminated surface than in the bulk, and about twice the 
bulk ferroelectric well depth. 

V. IN-PLANE POLARIZED PbTiOs (001) FILMS 

In this section we present our results for the c(2x2) 
PbTiOa (001) film with spontaneous polarization paral- 
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a) 



b) 




FIG. 4: Optimized atomic structure of the in-plane polarized 
c(2x2) (001) PbTi0 3 film with PbO surface termination, (a) 
Top view; only the first top two layers are shown, (b) Side 
view; only the top half of the 11-layer slab is shown. The 
oxygen atoms in the second layer are shown in a darker shade 
of gray. The in-plane Pb off-centering as well as the c(2x2) 
AFD O-rotation in the Ti02 subsurface layer, are clearly vis- 
ible. Pb-O bonds shorter than 2.6 A are shown in black. 

lei to the surface, P=P(100). In particular we are in- 
terested in studying the interaction between AFD and 
in-plane FE distortions in proximity of the (001) surface. 
We consider only the PbO termination, which is the most 
energetically favorable and, as we have shown in the pre- 
vious section, displays a strong surface AFD instability 
in absence of FE distortions. 

The optimized atomic structure of the c(2x2) in-plane 
polarized film is shown in Fig. Clearly visible are: (i) 
the in-plane Pb off-centering, present throughout the film 
and enhanced at the surface, and (ii) the large AFD Ti04 
rotation in the TiC>2 subsurface layer. Thus, contrary to 
the bulk where FE distortions suppress the AFD rota- 
tion, FE and AFD distortions coexist at the PbO surface. 

We now consider these coupled distortions in more 
detail, using a symmetry analysis. The c(2x2) super- 
cell can accommodate atomic distortions of the (lxl) 
ideal unit cell with wavevectors q=(000) and q=(110)7r/a 
(r and M point of the (lxl) surface BZ, respectively). 
Thus, the atomic displacements relative to the ideal per- 
ovskite structure can be decomposed into T- and M-point 
irreducible representations of the reference (lxl) unit 
cell: 8 a .i = $a i (where v labels the irreducible rep- 
resentations, i the atoms in the c(2x2) unit cell, and 
a=x, y, z). This decomposition allows us to separate the 
FE, AFE, and AFD contributions to the lattice distor- 
tion. We find that the in-plane atomic displacements can 
be decomposed into three irreps M 3+ , M 5 _, and T 5 (or 
E u ), which correspond to AFD, AFE, and FE distor- 
tions, respectively. The atomic displacement patterns 



for the M$ + , Ms_, and Ts_ irreducible representations 
are shown in Fig.[3^a-c). 

To quantify, layer- by-layer, the strength of the AFD, 
AFE, and FE distortions we proceed as follows. The 
layer-by-layer strength of the AFD contribution is de- 
scribed by the Ti04 rotation angle, 9, in each Ti02 layer 
(see Fig. El). For the FE (^=r 5 _) and AFE (u=M 5 -) 
contributions we define the average polar distortion along 
the a-direction, p 1 ^, in each (lxl) unit cell of each atomic 
layer (see Fig. EM); f° r Ti02 layers 

rt = «Ti)-MM±«M, (2) 

and for PbO layers 

_ [%(Pbi) + %(Pbn)] _ 

Pa — 2 °aV J )^ W 

where b v a are the atomic displacements in the a — x,y 
direction relative to the ideal (lxl) perovskite structure, 
for each irrep v. For a film with spontaneous polarization 
along x, the average FE polar distortion is along x, p^ E , 
and has the same sign in each (lxl) unit cell, while the 
average AFE polar distortion is along y, Py FE , and has 
opposite sign in each neighboring (lxl) unit cell. 

The computed values for the parameters describing the 
strength of FE, AFE, and AFD distortions in each atomic 
layer are summarized in Table IVll The FE distortion is 
enhanced at the surface and rapidly resumes the bulk 
value in the inner layers. The average polar distortion 
in the first PbO layer is 1.7 times larger than the cor- 
responding quantity in the bulk, and the enhancement 
is confined to the surface layer. The same behavior for 
the in-plane FE distortion was found at the (lxl) PbO- 
terminated surface^ where AFD distortion is not allowed 
by symmetry. There, the average in-plane polar distor- 
tion in the first surface layer was found to be 1.5 times 
larger than the corresponding quantity in the bulk. It 
should be noted that in that calculation, the in-plane 
lattice parameters of the slab were taken as the bulk 
tetragonal c and a parameters, and the comparison is 
to the tetragonal bulk polarization. In our calculation, 
we fixed the in-plane lattice parameters to match a sub- 
strate of bulk cubic PbTi03. The proper reference bulk 
limit in this case is the cubic lattice optimized allowing 
a spontaneous polarization along the (100) direction (i.e. 
allowing internal atomic displacements along (100)). The 
resulting "cubic polarized" structure (C-FE) is the bulk 
reference limit in Table IVTI 

The AFD oxygen-rotation is also enhanced at the sur- 
face, 02=10.8°, and rapidly vanishes in the deeper layers. 
Thus, the AFD distortion is confined to the surface oc- 
tahedron layer, and is comparable to the one found in 
absence of FE distortions. 

Finally, there is a small AFE distortion at the PbO 
surface layer that rapidly decreases in the inner layers 
where it is essentially zero. 
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TABLE VI: The average layer-by-layer ferroelectric, p x , and 
antiferroelectric, Py FE , polar distortions per (lxl) unit cell, 
in percentage of the bulk lattice parameter (at = 7.37a.u.), 
and the anti-ferrodistortive oxygen rotation, # AFD , in the 
Ti02 layers. The average spontaneous polar distortion and 
the AFD angle computed for the bulk "polarized cubic" struc- 
ture (C-FE) and the bulk "cubic AFD" structure (C-AFD) 
are given as a bulk reference. 



,AFE flAFD 



layer 


Ti0 2 


PbO 


TiO a 


PbO 


Ti0 2 


1 




11.9 




-2.4 




2 


4.1 




0.5 




10.8° 


3 




6.3 




0.4 




4 


3.7 




-0.06 




-0.1° 


5 




7.1 




-0.2 




6 


3.9 




0.08 




0.9° 


bulk (C-FE) 
bulk (C-AFD) 


4.1 
0.0 


7.1 

0.0 






0° 
3.3° 



The FE and AFD distortions have both the effect of 
shortening some of the Pb-0 bonds. At the surface, the 
compounded effect of the coexisting large FE and AFD 
distortions results in even shorter Pb-0 bonds, as short 
as 2.31 A (comparable to the shortest bondlength in PbO 
bulk). 

VI. CONCLUSIONS 

Non-polar cubic PbTiOs bulk has a weak AFD insta- 
bility, involving rotations of oxygen octahedra, which is 
suppressed by the FE distortion. In this work we have 
studied the effect of surface proximity upon the AFD 
instability at the (001) surfaces of PbTiOs. In partic- 
ular we studied the surface effect upon the instability's 
strength and its competition with in-plane ferroelectric- 
ity 

In the absence of FE distortions we found that the im- 
pact of the surface upon the strength of the AFD insta- 
bility is termination-dependent. At the PbO-termination 
the AFD deformation is strongly enhanced relative to the 
bulk; the top surface octahedra rotate by 11° about the 
surface normal. The enhanced AFD distortion is con- 
fined to the octahedron layer at the top of the surface. 
We estimated that the energy gain associated with the 
octahedron rotation is two orders of magnitude larger at 



the surface than in the bulk. At the Ti02-termination, 
the AFD deformation is essentially suppressed; the top 
surface octahedra (which are truncated) rotate by an an- 
gle of 0.1°. For both terminations the surface effect upon 
the AFD distortion is confined to the first octahedron- 
layer at the surface, and below the surface octahedra the 
AFD rotation recovers the small bulk-like value, « 3°, 
computed for non-polar cubic PbTi03. 

For polarized films (with in-plane FE distortion) we 
found that FE and AFD distortions coexist in the prox- 
imity of the PbO-terminated surface, leading to a struc- 
tural phase not found in bulk. The FE distortion is en- 
hanced at the surface layer as was previously found for 
the (lxl) films where the AFD distortions are not al- 
lowed by symmetry. Concurrently, the AFD distortion 
is enhanced at the surface as much as found in absence 
of FE distortions. Thus, AFD and in-plane FE distor- 
tions are both enhanced at the PbO surface and do not 
strongly interact. The enhanced AFD deformation does 
not suppress the FE distortion, nor is the AFD deforma- 
tion suppressed by the FE distortion. 

The surface effects are confined to the first perovskite 
unit cell at the surface (i.e. in the first three atomic 
layers), while bulk- like behavior is recovered below it. 

In summary, we have shown that in PbTi03 the 
strength of the AFD instability, consisting of rotation of 
the oxygen octahedra, and its interaction with ferroelec- 
tricity are affected by the proximity to the surface. At 
the PbO termination this leads to a phase, not observed 
in the bulk, where enhanced AFD and FE distortions co- 
exist. The enhancement of the AFD instability at the 
PbO termination is driven by the formation of shorter 
PbO bonds and is a consequence of the particular chem- 
istry of Pb; no enhanced AFD distortion is found, for 
example, at the surface of BaTi03. We suggest that the 
enhancement of the AFD instability at the PbO termi- 
nation should be common among Pb based perovskite 
oxides. 
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